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LONG-TERM  GOALS 

The  long-term  scientific  goal  of  my  research  is  to  better  understand  the  distribution  of  phytoplankton 
in  the  world’s  oceans  through  remote  sensing  their  influence  on  the  optical  properties  of  the  water.  An 
associated  goal  is  the  understanding  of  the  absorption  and  backscattering  properties  of  marine  particles 
in  terms  of  the  distributions  of  their  size,  shape,  and  composition. 

OBJECTIVES 

The  inherent  optical  properties  (IOPs)  of  marine  particles  are  most-often  modeled  as  homogeneous 
spheres  using  Mie  Theory.  Although  this  approach  has  been  fruitful,  the  next  logical  step  in  modeling 
marine  particles  is  to  abandon  the  normally-employed  spherical  approximation  and  use  more  realistic 
approximations  to  their  shape.  The  advent  of  computer  codes  capable  of  handling  more  complex 
shapes,  and  the  increased  computational  speeds  now  available,  suggest  that  particle  modeling 
employing  simple  non-spherical  shapes,  e.g.,  disks,  rods,  etc.,  could  become  routine.  For  example, 
Gordon  and  Du  (2001)  used  a  two-disk  model  to  try  to  reproduce  the  backscattering  by  coccoliths 
detached  from  E.  huxleyi  in  the  blue-green  region  of  the  spectrum.  They  found  that  the  gross 
morphology  of  the  particle  (e.g.,  modeling  the  coccolith  as  parallel  disks  as  opposed  to  a  single  disk  of 
the  same  diameter  and  volume)  was  paramount  in  determining  the  spectral  variation  of  backscattering. 
Later  Gordon  et  al.  (2009)  showed  that  the  resulting  spectral  variation  of  the  backscattering  cross 
section  for  the  two-disk  model  agreed  with  experiment,  and  that  additional  fine  structure  evident  in 
SEM  images  of  coccoliths  had  only  a  minor  influence  on  the  backscattering.  In  another  study,  Clavano 
et  al.  (2007)  looked  at  the  IOP  of  spheroids  (ellipses  of  revolution)  as  a  function  of  their  size  and 
aspect  ratio  (larger  major  axis  divided  by  smaller  major  axis)  using  refractive  indices  in  the  range 
found  for  marine  particles.  Their  results  suggest  that,  when  compared  to  equal  volume  spheres,  shape 
can  be  a  significant  factor  in  the  IOPs,  especially  backscattering. 

The  work  described  here  represents  the  culmination  of  my  research  into  computation  of  the  IOPs  of 
non-spherical  objects.  Of  particular  interest  to  me  was  the  scattering  and  absorption  of  particles 
showing  significant  deviation  from  spheres,  e.g.,  cylinders  with  large  aspect  ratios.  This  gross 
morphology  resembles  long  chain  phytoplankton.  My  objective  was  to  provide  answers  to  the 
following  questions.  How  do  their  scattering  and  absorption  properties  depend  on  the  length,  diameter 
and  index  of  refraction  of  the  cylinders,  or  more  significantly,  their  aspect  ratio  (AR  =  length  4- 
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diameter)?  Do  their  IOPs  reach  a  limit  in  which  they  become  simply  proportional  to  the  cylinder’s 
length,  so  their  length-normalized  IOPs  become  independent  of  aspect  ratio?  How  do  they  depend  on 
the  distribution  of  absorbing  pigments  within  the  cylinder?  How  do  they  depend  on  orientation?  How 
do  they  compare  to  the  scattering  and  absorption  by  spheres  containing  the  same  volume  of  material 
and  quantity  of  absorbing  pigment?  Is  the  spherical  approximation  effective  in  interpreting  their 
scattering  and  absorption  in  terms  of  an  equivalent  refractive  index?  My  goal  was  to  develop  an 
understanding  at  a  level  that  would  allow  the  incorporation  of  shape  infonnation  into  IOP  models  and 
in  the  analysis  of  IOP  data. 

APPROACH 

Of  the  several  computer  codes  that  allow  accurate  computation  of  the  scattering  and  absorption 
properties  of  individual  particles  with  smooth  nonspherical  shapes,  I  employed  the  discrete-dipole 
approximation  (DDA,  Draine,  1988;  Draine  and  Flatau,  1994),  as  I  had  used  it  earlier  in  my  study  of 
particles  similar  to  E.  huxleyi  coccoliths.  For  the  later  parts  of  the  study  I  used  an  improved  version  of 
the  code  (DDSCAT  7.0).  I  examined  randomly  oriented  finite  cylinders  with  length  up  to  25  pm  and 
diameters  up  to  1.5  pm.  At  this  size,  an  equal  volume  sphere  has  a  diameter  of  about  4  pm. 

WORK  COMPLETED 

I  have  completed  a  study  of  the  scattering  and  backscattering  of  non-absorbing  randomly-oriented 
cylinders  (diameter  of  0.5- 1.5  pm)  as  a  function  of  the  length  and  the  index  of  refraction  (1.02  to  1.20). 
The  goal  was  to  understand  at  what  aspect  ratio  the  backscattering  becomes  essentially  proportional  to 
the  length.  This  was  mostly  carried  out  using  the  DDSCAT  6.1  code.  For  the  study  of  absorbing 
particles,  I  have  modified  the  new  DDSCAT  7.0  code  to  include  cylinders  with  internal  structure,  i.e.,  a 
coated  cylinder.  This  structure  was  chosen  to  represent  a  long  chain  of  phytoplankton,  e.g.,  diatoms,  in 
which  the  absorbing  pigment,  if  any,  is  packaged  within  the  internal  cylinder. 

I  considered  absorbing  cylinders  ranging  in  diameter  from  0.5  to  1.5  pm  and  lengths  from  0.5  to  15 
pm.  Two  sets  of  cylinders  that  I  examined  were  homogeneous  and  had  refractive  indices  1.05-0.002/ 
and  1.05-0.010/.  Two  other  sets  were  coated  cylinders  in  which  the  inner  cylinder  had  a  diameter 
equal  to  half  the  cylinder  diameter  and  indices  of  1.05-0.008/  and  1.05-0.040/.  The  outer  cladding 
had  an  index  of  1.05  (no  absorption).  The  cylinder  with  homogeneous  index  of  1.05-0.010/  and  the 
coated  cylinder  with  core  index  1.05-0.040/  have  exactly  the  same  quantity  (mass  or  number  of 
molecules)  of  absorbing  pigment. 

The  quantities  that  I  computed  were  the  orientationally-averaged  extinction  cross  section,  <jc,  the 
scattering  and  backscattering  cross  sections,  cr/,  and  <Jbb,  the  absorption  cross  section,  <ja,  and  the 
scattering  phase  function  (differential  scattering  cross  section  divided  bycr/,).  The  extinction  and 
absorption  efficiencies,  Qc  and  Qa,  are  the  associated  cross  sections  divided  by  the  orientationally 
averaged  projected  area  of  the  particle  (the  area  of  the  particle’s  shadow  formed  by  the  incident  beam, 
ZD/2). 

The  final  report  of  this  work  was  submitted  to  ONR  April  15,  201 1. 
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RESULTS 


The  main  finding  of  the  study  was  that  for  AR  >  ~  3-5,  Qc  and  Qa,  become  independent  of  AR  and,  in 
fact,  close  to  those  computed  for  infinite  cylinders.  Similar  results  are  obtained  for  the  scattering  phase 
function  and  the  backscattering  probability  ( <Jbbl<Jb )• 

Figure  1  provides  the  computed  values  of  Qc  and  Qa  for  the  all  aspect  ratios  examined  with  absorbing 
cylinders  (m,  >  0)  and  in,-  =  1.05.  The  notation  in  the  legend  is  provided  in  the  figure  caption.  The 
upper  two  panels  are  for  all  aspect  ratios  (1/3  -  30)  and  the  lower  panels  for  AR  >  3.  The  lines 
correspond  to  exact  computations  for  infinite  cylinders  with  the  same  diameter  computed  with  the 
IPHASE  code  (Evans,  1988;  Fournier  and  Evans,  1996).  For  the  case  with  weaker  absorption,  the 
results  clearly  show  that  the  efficiencies  closely  follow  those  for  infinite  cylinders  as  long  as  AR  >  3, 
for  both  the  coated  and  the  homogeneous  cylinders.  For  the  case  with  stronger  absorption,  the  so- 
called  “package  effect,”  the  decrease  in  absorption  when  the  absorbing  molecules  are  not  unifonnly 
distributed  within  the  particle,  results  in  lower  absorption  efficiency. 

This  packaging  effect  displayed  more  clearly  in  Figure  2,  for  which  I  provide  the  ratio  of  the 
absorption  efficiencies  Qa  for  a  coated  cylinder  with  indices  minside  =1.05-0.040/  and  moutside 
1.05-0.000/  (^(Packaged))  to  that  of  a  homogeneous  cylinder  mmsuie  =  moutSide  =  1.05-0.010/ 
(^fl(Homo)).  Note  that  both  cylinders  contain  the  same  number  of  absorbing  molecules.  The  figure 
shows  that  the  effect  of  the  absorbing  pigment  packaging  is  greatest  in  the  blue  region  of  the  spectrum 
and  for  larger-diameter  cylinders.  The  maximum  decrease  in  Qa  due  to  the  packaging  is  about  25%  in 
this  case.  Although  the  symbols  do  not  differentiate  between  cylinder  lengths,  for  a  given  diameter  the 
packaging  effect  is  largest  in  the  shortest  cylinder  and  depends  very  little  on  the  length  once  AR 
exceeds  unity.  In  the  case  with  less  overall  absorption,  i.e.,  ///„„,*.  =1.05-0.008/  and  moutside  = 
1.05-0.000/  compared  to  that  of  a  homogeneous  cylinder  mmsuie  =  moutside  =  1.05-0.002/',  similar 
results  are  obtained  (not  shown);  however,  the  maximum  decrease  in  Qa  due  to  the  packaging  is  only 
about  10%. 

The  extinction  efficiencies  of  finite  cylinders  with  larger  values  of  p  are  compared  with  those  of 
infinite  cylinders  in  Figure  3.  In  this  case  the  refractive  index  is  1.20,  and  since  there  is  no  absorption, 
Qb=  Qc  •  The  figures  clearly  show  that  for  AR  >  3,  the  scattering  efficiency  is  again  close  to  that  of  an 
infinite  cylinder,  with  all  cases  except  two  differing  from  an  infinite  cylinder  by  less  than  ±10%  (RMS 
difference  ~  5%). 

I  have  shown  that  the  extinction  and  absorption  efficiencies  of  micrometer-sized  cylindrical  particles 
depend  little  on  the  aspect  ratios  as  long  as  AR  >3.  Is  this  the  case  for  the  phase  function  and 
backscattering  probability?  Figure  4  provides  an  example  of  the  weak  dependence  of  the  scattering 
phase  function  (and  degree  of  linear  polarization)  on  aspect  ratio.  Virtually  the  only  differences  in  the 
phase  function  between  AR  =  20  and  AR  =  3  are  the  enhanced  scattering  near  zero  degrees  and  the 
deeper  minima  near  30°,  50°,  80°,  and  130°  for  AR  =  20.  This  weak  dependence  on  aspect  ratios  as 
long  as  AR  >  5  is  also  displayed  by  the  backscattering  probability  as  shown  in  Figure  5  for  a  wide 
range  of  refractive  indices  and  particle  diameters. 
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IMPACT/APPLICATIONS 


The  fact  that  (normalized)  IOPs  by  cylinders  with  aspect  ratios  of  3-5  become  independent  of  the 
aspect  ratio  can  be  used  to  simplify  computation  of  such  for  analysis  of  particle  scattering.  For 
example,  computation  of  the  scattering  efficiencies  for  infinite  cylinders  (which  can  be  carried  out  very 
efficiently)  can  be  used  to  estimate  the  same  for  cylinders  with  AR  >  3  with  an  error  <  5-10%.  These 
efficiencies  can  then  be  used  to  compute  the  cross  sections  for  the  cylinder  length  of  interest.  Similarly, 
the  backscattering  probability  for  long  cylinders  with  diameters  ~  1-3  times  the  wavelength  (which 
cannot  be  computed  with  most  light  scattering  codes)  can  be  deduced  (within  5-10%)  from  those  of 
finite  cylinders  with  AR  ~  3-5.  These  observations  can  simplify  the  analysis  of  light  scattering  and 
absorption  by  high-aspect  ratio  particles. 

RELATED  PROJECTS 

None. 
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Figure  1.  Extinction  and  absorption  efficiencies  computed  for  randomly  orientated,  homogeneous 
or  coated  cylindrically  shaped  particles,  given  that  their  diameter  and  aspect  ratio  are  known  so  that 
their  orientationally-averaged projected  area  is  DL/2.  Here,  p  =  2a  (m,  -l)  and  p'=  4antj,  where 
mr  -im„  is  the  refractive  index  of  the  particle  relative  to  water,  and  a=  ttD/X,  with  D  the  cylinder’s 
(outer)  diameter  and  X  the  wavelength  of  light  in  the  water.  Solid  lines  are  the  exact  computations 
for  randomly  oriented,  infinite  cylinders.  The  notation  “m=1.05-040i_1.05-000i”  indicates  that  the 
refractive  index  of  the  core  is  1.05  -  0.040i,  and  the  refractive  index  of  the  coating  is  1.05  -  O.OOOi, 
etc.  In  the  case  of  coated  cylinders,  p  and  p'  are  computed  using  nij  of  the  associated  homogeneous 
particle.  Top:  all  aspect  ratios  (1/3  -  30).  Bottom:  all  aspect  ratios  >3.  The  graphs  clearly  show  that 
the  efficiencies  closely  follow  those  for  infinite  cylinders  as  long  as  AR  >3,  for  both  the  coated  and 
the  homogeneous  cylinders.  For  the  case  with  stronger  absorption,  the  so-called  “ package  effect,  ” 
the  decrease  in  absorption  when  the  absorbing  molecules  are  not  uniformly  distributed  within  the 

particle,  results  in  lower  absorption  efficiency. 
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Figure  2.  This  figure  provides  the  ratio  of  absorption  efficiencies  (coated  to  homogeneous)  of  the 
more  strongly  absorbing  cylinders  as  a  function  of  wavelength.  The  diameter  of  the  cylinder  (in  pm) 
is  specified  in  the  legend.  For  each  diameter,  the  symbols  refer  to  cylinder  lengths  ranging  from  0.5 
to  15  pm.  The  figure  shows  that  the  effect  of  the  absorbing  pigment  packaging  is  greatest  in  the 
blue  region  of  the  spectrum  and  for  larger-diameter  cylinders.  Although  the  symbols  do  not 
differentiate  between  cylinder  lengths,  for  a  given  diameter  the  packaging  effect  is  smallest  in  the 
shortest  cylinder  and  depends  very  little  on  the  length  once  the  aspect  ratio 
(length/diameter)  exceeds  unity. 
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Figure  3.  Qc  computed  for  non-absorbing  cylinders  (m  =  1.20  -  O.OOOi)  with  diameters  (D)  ranging 
from  0.5  pm  to  2.0  pm.  Left:  0.25  <AR  <30  (points  colored  in  red  are  AR  =  2).  Right:  3  <AR  < 
30.  The  solid  curve  is  the  extinction  efficiency  (for  a  unit  length)  of  randomly-oriented  infinite 
cylinders.  As  in  Figure  1,  p  =  2a  (m  -1)  with  a  =  zrD/A.  The  figures  clearly  show  that  for  AR  >3, 
the  scattering  efficiency  is  again  close  to  that  of  an  infinite  cylinder,  with  all  cases  except  two 
differing  from  an  infinite  cylinder  by  less  than  ±10%  (RMS  difference  ~  5%). 
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Figure  4.  Orientationally-averaged  scattering  phase  functions  (left)  and  degree  of polarization 
(right)  at  600  nm  (vacuum)  for  homogeneous  cylinders  with  a  diameter  of  1  pm  and  length  of  3  pm 
(DlxL3)  and  20  mm  (DlxL20).  The  refractive  index  is  1.05  -  0.0021  Virtually  the  only  differences 
in  the  phase  function  between  AR  =  20  and  AR  =  3  are  the  enhanced  scattering  near  zero  degrees 
and  the  deeper  minima  near  30  °,  50  °,  80  and  130  °for  AR  =  20. 
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Figure  5.  Examples  of  the  variation  of  the  backscattering  probability  with  aspect  ratio  for  cylinder 
diameters  between  0.5  and  1.5  pm  and  refractive  indices  ranging  from  1.02  to  1.20.  The  black 
curves  are  for  a  vacuum  wavelength  of 400  nm  and  the  red  curves  for  700  nm.  The  graphs  show 
that  there  is  very  little  variation  in  the  backscattering  probability  with  aspect  ratio  for  AR  >  5. 
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